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We construct phantom energy models with the equation of state parameter w which is less than 
— 1, w < —1, but finite-time future singularity does not occur. Such models can be divided into two 
classes: (i) energy density increases with time ( "phantom energy" without "Big Rip" singularity) and 
(ii) energy density tends to constant value with time ( "cosmological constant" with asymptotically 
■ de Sitter evolution). The disintegration of bound structure is confirmed in Little Rip cosmology. 

7-H ' Surprisingly, we find that such disintegration (on example of Sun-Earth system) may occur even in 

asymptotically de Sitter phantom universe consistent with observational data. We also demonstrate 
that non-singular phantom models admit wormhole solutions as well as possibility of Big Trip via 
^ " wormholes. 

<D '• 

t^H , Introduction. The discovery of accelerated expansion of the universe [J [|| led to a number of new ideas/solutions 
in cosmology. Recent observations of supernovae are consistent with the universe made up 71.3% of dark energy 

t-H [ and only 27.4% of a combination of dark matter and baryonic matters Q- Dark energy proposed to explain the 
cosmic acceleration should have the strong negative pressure (acting repulsively) in order to explain the observed 
\ accelerating expansion of the universe (for recent reviews, see |4j-l9|)- The equation of state parameter wd for dark 
energy is negative: 

kJj w D =p D /p D <0, (1) 

where pu is the dark energy density and p^> is the pressure. We omit subscript D further for simplicity. 

According to the latest cosmological data available, the uncertainties are still too large to discriminate among the 
three cases w < — 1, w = — 1, and w > — 1: w = — 1.04t^io US EH- If w < —1, the violation of all four energy 
conditions occurs. The corresponding phantom field, which is unstable as quantum field theory [l2T | but could be stable 
in classical cosmology, may be naturally described by a scalar field with the negative kinetic term. Such Lagrangians 
appear in some models of supergravity [l3T |. in the gravity theories with higher derivatives flij and in string field 
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The additional interest to the models with the phantom fields is caused by their prediction of a so-called Big Rip 
^ singularity [l6l - [2l| . Theoretically, the scale factor of the universe becomes infinite at a finite time in the future which 
1 _— I . was dubbed Big Rip singularity. There were proposed several scenarios to cure the Big Rip singularity: (i) To consider 
phantom acceleration as transient phenomenon. This is possible for a number of scalar potentials, (ii) To account 
for quantum effects which may delay/stop the singularity occurrence [22j . (iii) To modify the gravitation itself in 
such a way that it appears to be observationally-friendly from one side but it cures singularity (for review, see f23j). 
(iv) To couple dark energy with dark matter in the special way [24| or to use special (artificial) form for dark energy 
equation of state [25j . Note that for quintessence dark energy, other (milder) finite-time singularities may occur. The 
corresponding classification of such quintessence-related finite-time singularities is given in Ref. 34]. For instance, 
type II (sudden) singularity [H[ or type III singularity [34| occurs with finite scale factor but infinite energy and/or 
pressure. Such quintessence-related finite-time singularities occur for instance, for the models (2fj|-l30| and were also 
called the "big freeze" singularity [U H2| • 

The closer examination shows that the condition w < — 1 is not sufficient for a singularity occurrence. First of all, a 
transient phantom cosmology is quite possible. Moreover, one can construct such models that w asymptotically tends 
to —1 and the energy density increases with time or remains constant but there is no finite-time future singularity 
[33T - r38j . Of course, most evident case is when Hubble rate tends to constant (cosmological constant or asymptotically 
de Sitter space), which may also correspond to the pseudo-rip |39j. Very interesting situation is related with Little 
Rip cosmology |38j where Hubble rate tends to infinity in the infinite future (for further investigation, see [39ll40l.[5l|). 
The key point is that if w approaches —1 sufficiently rapidly, then it is possible to have a model in which the time 
required for singularity is infinite, i.e., the singularity effectively does not occur. Nevertheless, it may be shown that 
even in this case the disintegration of bound structures takes place in the way similar to Big Rip. 

The aim of this article is to develop the method of constructing the phantom models without finite-time singularity. 
In Sec. II, a general approach to this problem is developed. The examples of singular dark energy models are given 
there. Sec. Ill is devoted to the construction of the scalar Little Rip dark energy models. In Sec. IV, the transient 
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phantom era which ends up at asymptotically de Sitter universe is investigated. The corresponding non-singular 
scalar phantom models are constructed. It is demonstrated that such models are compatible with latest data from 
Supernova Cosmology Project. We show that the dissolution of bound structures is possible in such asymptotically de 
Sitter universe for special choice of theory parameters. This gives the observationally-consistent pseudo-rip cosmology 
scenario. The influence of possible interaction between phantom energy and dark matter on non-singular cosmological 
evolution is investigated in Sec. V. In Sec. VI, the possibilities of wormhole solutions and so-called "Big Trip" 
wormhole scenario in constructed cosmological models are considered. Some summary and outlook are given in the 
Discussion section. 



Scalar dark energy models with future singularity. We start from the FRW equation and the conservation law 
for spatially flat universe 

dV = f, P = -3(£)(,+p). (2) 

where p and p are the total energy density and pressure, a is the scale factor, ' = d/dt, and we use the natural system 
units in which 8ttG = c = 1. 

We will examine the future evolution of our universe from the point at which the pressure and the density are 
dominated by the dark energy. For the pressure of dark energy, one can choose the general expression 

p = -p- f(p) , (3) 

where f(p) is a function of the energy density. The case f(p)>0 corresponds to w < — 1. From ([2]), one can obtain 
the following relation between the time coordinate t and f(p): 

2 f* dx 

t = TzLW)' x = ^ p - (4) 

If the integral Q converges at p — > oo, we have a singularity: energy density becomes infinite at a finite future t = tf 
in the future. The expression for scale factor 

(2 r xdx\ 

a = a exp - / — — , (5) 



3 J Xo f{x) 

indicates that there are two possibilities: 

(i) scale factor diverges at a finite time ("Big Rip"). 

(ii) scale factor reaches finite value and singularity (p — > oo) occurs. It is type III singularity in the notations of 
Ref. [13. 

The simple choice of / corresponding to first scenario is 

f{x) = Px a , l<a<2, (6) 

where j3 is a positive constant. The scale factor can be written as 

Aexp{(B - Ct)-~<} , a 7 = fcr m 

[D-Et)- S , a = 2, 5 = 2/3/3, 1 ' 

where A, B, C, D, and E are positive constants. The case a — 2 corresponds to a simplest model of phantom energy 
with parameter w = — 1 — j3 = const. 

If a > 2, one has the second possibility: the energy density grows so rapidly with time that scale factor does not 
reach the infinite value. 

Equivalent description in terms of scalar theory can be derived using the equations: 

P=-<j> 2 /2 + V(<P), p=-4> 2 /2-V{^) (8) 
where is a scalar field with potential V (</>). For the scalar field and the potential, one can derive following expressions: 

2 f x dx 

V(x) =x 2 + f(x)/2. (10) 
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Combining Eqs. © and (fT0|) . we find the potential as a function of scalar field. For simplicity we choose the sign "+" 
in Eq. ©, hereafter. The choice ^ yields for potential (if a^2) 

V(4>) =F^(4>-<t> )^7* +^1^72(^-^)1^72 , F = y/zp{\-a/2)/2. (11) 
For a = 2, as expected, we have an exponential potential: 

V{<t>) = (1 + 0/2)4 exp{ (0 - O )} • (12) 

The key difference between (i) and (ii) for this model is that the potential of scalar field has a pole in a case of 
(ii)-singularity. Note that another type of singularity occurs if f(x) — > oo at x — Xf < oo, i.e., the pressure of dark 
energy becomes infinite at finite energy density. The second derivative of scale factor diverges. 

It is interesting however to investigate the phantom energy models without finite-time future singularities with the 
help of similar technique. 



Scalar Little Rip cosmology. Let us now consider the models which provide an evolution for the universe inter- 
mediate between de Sitter evolution and phantom era with the Big Rip. These models were described in detail in 
[HI, HH, [5l| . The energy densi ty g rows with time but not rapidly enough for the occurrence of the Big Rip singularity. 
According to terminology in |38| we have a so-called "Little Rip" : eventually a dissolution of bound structures at 
some point in the future occurs. 

The method of constructing such models is very simple (compare with [5l[). In principle one can choose arbitrary 
monotonic function g(x) well defined in domain x > xq and satisfying condition g{x) — > oo at t — ¥ oo. Then one can 
assume that 



From Eqs. (j3| and (JSJ we have 



/(*) = ■ (13) 
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t = -^(g(x) - g(x )) , 

a = a Q expl^((xg{x)-x g(xo))~'^ J g(x)dxj . (14) 
For the scalar field <p and potential V(<p), the following expressions can be written 

0(x)=0 o + ^^ VTixjdx, V (x)=x 2 + ^-j. (15) 
The exponential growth of density with time (g(x) = ln(x)//3, (3 > 0) corresponds to the scalar potential 

V{4>) = ^ + ^<f> 2 , 4>= ^ exp {^-t\ , (16) 



(for simplicity, we put 4>q = A^xq/SP). The scale factor a grows with time in accordance with the double exponential 
law 

(2xo ( ( 3 1/2 ?t 



a = a exp j — (^exp j-ljj. (17) 

As shown in (38[, this model can be matched with the latest data from Supernova Cosmological Project. The best- fit 
value for (3 is 3.46 x 10 -3 Gyr -1 . Thus, we presented the example of the scalar Little Rip cosmology where the future 
singularity does not effectively occurs. 

However, it should be noted that the Little Rip produces the disintegration of bound structures just as in the case 
of Big Rip. The condition of disintegration can be derived in the following. The acceleration of the universe leads to 
an inertial force on a mass m as seen by a gravitational source separated by a comoving distance / 

F ln = ml- = ml^—{2p(a) +p'(a)a). (18) 

(X o 



4 



The structure disintegrates when the inertial force (fl"8|) . dominated by dark energy, becomes equal to the force 
bounding the structure. It is convenient to define dimensionless parameter [39j 

F ia = 2p(a) + //(a)Q . (19) 
Po 

(po is a dark energy density at the present time) . The simple calculations allow to derive the following expression for 
Fi n as function of time 

F in = 2 eMZ 1/2 Pt) + jr 2 exp {^f 1 ) ■ (20) 

The system of Sun and Earth, for example, disintegrates when F ln reaches ~ 10 23 . Therefore, the time required for 
this event is around 8.5 x 10 3 Gyr. 

Phantom models with asymptotically de Sitter evolution. Another interesting class of models arises if integral 
in Eq. Q diverges at some finite x — Xf < oo. The time required for energy density to reach p — x 2 is infinite, i.e., 
the expansion of the universe asymptotically approaches the exponential regime, which corresponds to the pseudo-rip 
in (3i| . The energy density tends to the constant value ( "cosmological constant" ) although the parameter w is always 
less than —1. 

For example, let us assume that 

f{x) = A(l - x/x f ) a , (21) 

where A and a are positive constants and we assume a > 1. In this case, the integral ((4]) diverges at x — Xf . For the 
case that q / 1,2, algebraic calculations allows us to get the following representation for scale factor: 

a(t) = a exp( Xl t/V3) exp(g a (t)) , g a (t) = - ~ — + ( X " ^) ) ■ ( 22 ) 

For 1 < a < 2 one can easily see that g(t) — > for t — > oo. If a > 2 \g(t)\ <C Xft/y/3 at t — > oo. Therefore the 
dependence (|2ip asymptotically tends to de Sitter solution with vacuum energy density p\ = xf . 
When t oo, the value of scalar field 

2xt 1 \( x y- a/2 [V3(l-a)At ' - ^-"^ 



3A 1 — a/2 I \ Xf ) \ Xf \ Xf 



I"?) ! >> ( 23 ) 



tends to constant for 1 < a < 2 and to — > ±oo for a > 2. As described in the second section, the scalar potential 
may be found as 



v(4>) = x\ i - ' ' 1 x ° 



'- a/2 V3A(l-a/2)(cf>-4> ) S 



1^72 



Xf J 2xf I 

a\l x v- a/2 V3A(i- a /2)( ( j ) -M \ T ^ r2 ^ (24) 



The cases a = 1, 2 are more interesting. For a = 1, the scale factor behaves as 

2x\ ( x 
3A 1 xf 



a(t) = a exp(xft/V3)ex V { gi (t)) , 9l (t) = ^l(l-^ \ (exp(-V3 At /2xf) - 1) , (25) 



and the scalar field 

1/2 



= <Po + ^jL ( 1 - ^ ) (1 - exp(- V3At/4x f )) , (26) 
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tends asymptotically to maximum (if xj > A/4) or minimum (if xj < A /4) of corresponding potential 

v{4>) = x\ + \(d> - <n 4 - m 2 (0 - <n 2 , r = O + l§= (1 - x - x ' ' 



3A \ x f 



9 A 2 2 _ 3 A 2 fix 2 



16xf 8x 2 \ A 

The choice a = 2 leads to the exponential potential 



* = 773 = "53 I " 1 ) ■ (27) 



a(t) = a Q cxp{x i t/V3)exp{g2{t)) , 3 2 W = -^j-lnl ^ i + 1 



(28) 



^Infl + ^V (29) 



/3A 
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V = xj- 2x f exp(-S) + (1 + ^2 J exp(~25) , S* = ^3^(0 - o )/2x f . (30) 

The appearance of exponential potentials may indicate to some connection with string theory. Eqs. (|29p and (|30[) 
show that in the infinite future, the scalar field also goes to infinity and the scalar field climbs up the potential to 
constant V — > x\. 

The above example is a good theoretical illustration of the new phantom energy models mimicking vacuum energy. 
The dark energy with such a behavior can be realized if the function f{x) is equal to zero at x = Xf and the integral (U]) 
diverges at x = Xf . In the vicinity of this point, the arbitrary function / satisfying these conditions can be expanded 
as 

f(x) = (x - Xi ) a + 0((x - x { ) a ) , a>l. 

Thus for t -> CO, one can expect that Eqs. (J23J, (J2BJI, and $29$ for <p(t) and (JMP, dHJ), and for V(<f>) will 
be satisfied. The case a = 1 corresponds to the most rapid growth of the energy density with time in which any 
singularity does not occur. 

Let us consider another model: 



/(*) = W^-j. (31) 

In the vicinity x = X{ , we find f(x) w ^ir(l — x/xf) 2 and we can conclude that the potential of scalar field looks like 
that in (|3"0"|) at i —> oo. Indeed, this is correct. For such a model, one gets 

V(</>) = x 2 [l arctanexp(S) 1 + , S=^(0-0 O ). (32) 

L 7T J (l + exp(6)) J Xf 

The interesting question is: could such models in principle describe latest supernova data from the Supernova Cos- 
mology Project? The analysis shows that it is possible. Moreover the construction of such models is trivial. For 
example let us choose 

f(x) = (3x^(1 -( X / Xi f/ 2 ), (33) 

and assume that dark energy density varies from to pi — x\ . Let t = to and p — po ai current universe. Eq. ([5]) 
allows to write the following relation between the dark energy density p and the redshift z = ao/a — 1: 

p(z)^p i (l-(l + zy'(l-A)f 3 , 7 = 3/3 Pf - 3/4 /2, A=^V /4 . (34) 
The equation of state parameter wq is 

2 7 1 - A 
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FIG. 1: The difference between modulus extracted from the toy dark energy model (|33[) and from the standard ACDM-model 



For the scale factor, we have the parametric expression 

f'o 



a(v) 
t(v) 



(1 - w 3/2)2/3 7 ' 

2V3 A, , , , 1 2v + l 7T 1, „ 

to + —[- 6 Hv +v + i) + ^ - - H i- V) 



2/3 7 



ao exp 



Xf(t - t Q ) 



For the dark and baryonic matter densities, 

Prn = Pmo(l + z) 3 . 

Therefore, the dependence of luminosity distance Dl from redshift z for this model is given by 



Dr 



(i + z) / (n m (i + zf + OdHz)) 



-1/2 



dz . 



h(z) = A~ 4 / 3 (l — (1 + z) 7 (l — A)) 



4/3 



(36) 
(37) 



The parameter v varies from (at t = to) to 1 (at t —> oo). The last term in Eq. (|3"T|) dominates at u — > 1. In this 
case, with a good accuracy one can write 



(38) 



(39) 



(40) 



In Eq. (|40p , fi m and J7d express the fractions of matters and dark energy in the total energy budget correspondingly. 
For the standard ACDM cosmology, as well-known, we find 



D[ C = + z ) f ( Q m(l + z f + n A ) 1/2 dz . 



(41) 



If 7 and A are very close to zero and 1 correspondingly (/3 is small and p — ¥ pi), then h(z) is close to 1. Therefore 
the functions (j40|) and (|4T|) are essentially indistinguishable (especially in the observable range < z < 1.5). The 
equation of state parameter is —1. Another choice of parameters 

A = 0.5, 7 = 0.075 

corresponds to current equation of state parameter wq = —1.05. The difference <5/i = b\g(D / D sc ) (/i is a distance 
modulus) for < z < 1.5 is depicted on Fig. [T] and does not exceed 0.016 (for fi TO = 0.28, ft a = = 0.72). Taking 
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FIG. 2: The time dependence of w for the model (|33 
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FIG. 3: The time dependence of the dark energy density for the model 
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into account that the errors in definition of SNe modulus are ~ 0.075 -j- 0.5, we conclude that our model fits these 
data with excellent precision. 

The equation of state parameter w and phantom energy density very slowly increase with time (Fig. [21 Fig. [3|). 

Thus, we presented non-singular phantom dark energy which evolves to asymptotically dc Sitter space and satisfies 
the observational bounds. 

Let us compare the model (|33|) with observational data in more detail. One defines deceleration parameter qo, jerk 
parameter jo as follows 



1 d 2 a 



qo = 



Jo 



aH 2 dt 2 
f 1 d 3 a 



1 



ld(H 2 ) h2 



t=t 



t=o 



H 2 1 2 dN 
1 d 2 (H 2 ) 3 d(H 2 



N=0 



2H 2 dN 2 



2H 2 dN 



(42) 



AT=0 
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Here N is the e-foldings defined by 

N = -hx(l + z). (43) 
For the current universe t = to, we have N — 0. Since 



one gets 



H 2 = \ (p + Pm ) = I (pf (1 - (1 - A)) 4/3 + Pm0 c- 3N ) , (44) 

Jo = Q7 2 (A- 1 - l) 2 + (-^ 7 2 + 2 7 ) (A- 1 - 1)) n A + 1 . (45) 

Since A = 0.5, 7 = 0.075, Vt m = 0.28, and A = 0.72, we find 

?0 = -0.54, Jo = 1.11. (46) 

In case of ACDM model, which corresponds to A = 0, we have qo = —0.58 and jo = 1, which is not so different from 
the values (|46]) . 

Just as it was done in previous section, we can estimate the possibility of the disintegration of bound system in our 
model (|3"3"|) . It is easy to derive the following parametric representation for the dimensionless inertial force F{ n in (flU|) 
as follows, 

F in (u) = 2A~ 4 / 3 (1 - (1 + U )?(l - A)) 1 ^ |i + g 7 _ ^ (1 + u) 7(i _ A )| , (47) 
1 f° 

t-t = — du{i + u)- 1 (n m (i + uf + n B h(u)V 1/2 . (48) 

The variable u = ao/a — 1 varies from (at present time) to —1 (when t — > 00). The function h(u) coincides with 
h(z) in Eq. (|40|) by changing z — !• u. The inertial force asymptotically tends to 2A~ 4 / 3 . Therefore, the disintegration 
of the system of Sun and Earth can occur if A < A m j n = 10~ 17 . The parameter 7 at A = A m i„ may vary from 
(when wo — —1) to 2 x 10~ 18 (when Wq = —1.14 ~ the lower bound from observations). The analysis shows that for 
any A < A m i„ and < 7 < 10^ — 1)A, our model also describes the SNe data. Therefore, the current observations 
in principle do not contradict with the possible disintegration of bound system in the models with an asymptotically 
de Sitter expansion. Note that such models were dubbed pseudo-rip in Ref. [39j | . 

One remark is in order. The asymptotically de Sitter expansion of the universe can occur not only in the models 
with the phantom energy but also for the quintessence dark energy (— 1 < w < 0). For the quintessence, the function 
f(x) in is negative and the energy density decreases with time. Eqs. (JU), ([5]), and (fit)]) do not change, for Eq. 
it is sufficient to replace f(x) — > —f{x). One can consider our model (|2Tj) and assume that a = 1 and xq > Xf. The 
energy density asymptotically tends to p — > x\ but the equation of state parameter to to > -1. Therefore one can 
conclude that the point x — Xf is attractor in the case f{x) = A(l — x/x{). 



Coupled phantom models. Realistic cosmological scenarios should take into account that the dark energy does 
not have a single component of universe energy. It is quite appealing to include the possible interaction between dark 
energy and dark matter (for recent discussion, see [4ll - l43l ]). We shall see that the addition of such a coupling may lead 
to interesting effects in the non-singular phantom cosmology. The cosmological solutions with interacting phantom 
energy exhibit much richer behavior than those considered above. 

It is customary to assume that phantom dark energy and dark matter interact through a coupling term Q as 

p m + 3H Pm = -Q , (49) 
p u + 3H{p D + PD ) = Q. (50) 

Let us introduce variables x and y 

PD=X 2 , p m =y 2 - 

For convenience, we assume that 



Q = ZHq{x,y) , 



(51) 
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where q(x, y) is an arbitrary function of variables x and y. Then the Hubble rate H is given by 

H=(^)"\ (52) 



and Eqs. (|49|) and (|50f can be rewritten as 

31/2(^2+ 2)1/2 

(2/ 2 +<?(z,y)) , (53) 



22/ 

31/2(^2+^2)1/2 

2x 



(/(*) +«(*,!/)) ■ (54) 



We can obtain critical points of the system (|53"j) . (f5~4"]) satisfying following conditions, if = and jff = 0. Note that 
these critical points must satisfy the requirement xt > 0, yt > 0. Unfortunately, there are no well-motivated direct 
observational or theoretical bounds on the type of interaction. Usually, simple interactions q ~ p m , q ~ po, and 
9 ~ (pd + Pm) are extensively investigated. 

For simplicity, we restrict ourselves to the cases q(x,y) = g(x), i.e., the intensity of interaction depends only from 
the phantom energy density and the rate of the Hubble rate H . Let us consider two examples, 

1. The simplest phantom energy model with f(x) = fix 2 and w = —1 — 0. Without interaction, one gets the Big Rip 
singularity. One can assume that 

q(x, y) — —ax 4 , a > , 

i.e., interaction leads to transformation of phantom energy into dark matter. On the other hand, the density of 
dark energy grows with time. Eventually the dynamical equilibrium between these processes is established. The 
equilibrium phantom energy density is 

en 

The density of dark matter otherwise decreases as universe expands but the dark matter born out from the phantom 
energy. The density of the dark matter asymptotically tends to 

pZ = — ■ 



a 



Then the Hubble rate tends to 



3 a 

i.e., we have expansion according to the de Sitter law with "effective" cosmological constant 

off £ 



a 

In other words, coupling with dark matter changes the qualitative behavior of phantom dark energy making it to be 
non-singular, in the way similar to the one described in Ref. [24| . 

2. The "switch" from the Little Rip to the asymptotically de Sitter expansion. As an example, let us choose 

f(x) = fix, q(x) = —ax 2 , a > . 

The interpretation of interaction is the same as in the previous model. Instead of the Little Rip expansion, the 
quasi-de Sitter expansion occurs with the effective cosmological constant: 

A eff = 4(a + 1)- 
or 

Thus, we demonstrated that coupling of phantom dark energy with dark matter may help in the transition from 
singular accelerating expansion to the non-singular one. Note that stability of such cosmologies may be investigated 
in the same way as in Ref. [44j where it is shown that the Little Rip often may be stable if we compare it with de 
Sitter universe. 
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"Big Trip" in phantom cosmology without singularities. One of the very interesting points in phantom 
cosmology is connected with wormholes. First of all, the existence of the static, spherically symmetric wormholc 
solutions of the gravitational field equations in the absence of ghost (or phantom) degrees of freedom is impossible, 
as it was shown in 14 511 . We discuss here wormholes application in phantom cosmology under discussion. 

It was shown in [4(1 |4?J that as one goes towards the Big Rip, there would occur the process of the fast wormhole 
swelling taking the size of the wormhole throat to infinity during the finite time. The reason of such a striking behavior 
is the phantom energy accretion onto wormhole. This accretion induces an increase of the wormhole throat radius 
so quick that the wormhole would engulf the entire universe before this reached the Big Rip. Such a result has been 
dubbed as "Big Trip" and was later criticized in Ref. [48| . The rejoinder was contained in [4{|. The issue remains 
open up to now and it is not proposed to deal at length with this discussion. All we need here is to consider the formal 
possibility of the Big Trip in phantom models without Big Rip. As we shall see, the Big Trip in phantom cosmology 
is a common occurrence, that often happens even if we consider models with asymptotic dc Sitter evolution. 

The equation which describes evolution of the throat radius of a Morris-Thorne wormhole b = b(t) due to dark 
energy accretion has the form 



-2ttW(1 + w)p, 



(55) 



where D is positive dimensionless constant and p is the energy density of the dark energy fluid. This equation is 
immediate consequence of the result for the dynamics of the mass of a black hole due to fluid accretion [50] ■ In this 
section we consider four examples of phantom models without finite-time singularity. 

1. The first and second model was obtained in (38ll5l|. Let put the Hubble rate as H — Hoe xt with positive constants 
Hq and A. As a result we have the scale factor, the density and the pressure in the form 



a(t) = ao exp 



H, 



(e At - 1) 



p{t) = ff V At , p(t) 



^e At [2A + 3ff e A <], 



(56) 



so 



w(t) = -1- 



2A 
3#o 



-At 



and at t oo one get the asymptotic de Sitter universe. 
Substituting (|5"6")) into (|55|) . one obtains 



Thus we have the Big Trip at 



b(t) 



4w 2 DH (e xt BT - e M ) 



(57) 



t = t r 



1 
A 



log b + 



4ir 2 DH Q 



where bo is the throat radius of the wormhole at t = 0: b = b(t). 

2. Let H = Ho — H\e~ xt . This is the most clear example of universe with de Sitter asymptotic behavior, if Hq > 
and A > 0. After integration of FRW equations we have 



J(t)=-1- 



2A# 1 c 



-A/ 



2(#o-#ie- A ') 2 ' 



and w(t) — > —1 at t — > +oo. It is interesting to note that this solution contains so called w-singularity [52j, to be 
more precise, some generalization of w-singularity which was obtained in [53j |. In fact, at 



t ty 



1 . H 
— log — 
A Hi 



w(t w ) — —oo although p(t w ) — and p(t w ) = —2\Hq/3 7^ 00. Sure, the w-singularity occurs if both Ho and Hi are 
positive constants. If Hi < then w(t) is always finite. 
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After calculations we have the throat radius of the wormhole in the form: 

3-4^ (1-e-*) - (58) 

For negative values of H± we have Big Trip for any values of positive parameters bo and D but have no w-singularity. 
If Hi > then Big Trip takes place only for wormholes with 

•} 

b > 



4tt 2 DHi ' 
at 

t = t BT = -\\og ( 1 - 



For positive Hi we have w-singularity and 

b(t w ) = 

If 



A b V ^ 2 DHibo 



36 

3 + 4 7 r 2 J D6 ( J ffo- J ffi) 



fli > Ho 



4tt 2 6 ' 



then Big Trip takes place after w-singularity, otherwise - before. In any case, it is clear that w-singularity does not 
influence the evolution of wormhole at all, since it is not real physical singularity. 

3. Now let consider the solution (0 with a^2. Since we are interested in solutions without Big Rip one should take 
7 = — g < 0; this is possible if a < 1 (g > 1) and a > 2 (0 < g < 1). We have phantom model w < —1 for the a < 1, 
and w > —1 for the a > 2. It is interesting to note that if g > 1 (a < 1, w < —1) then for t — >• i s = B/C we have 
p -> 0. If 1 < g < 2 (a < 0, w < -1) then at t -> t s p — !• — oo, p->0. This is type II singularity [34]. If g > 2 
(0 < a < 1, io < —1) then p(t s ) = p(t s ) = 0, w(t s ) = —oo. This is exact w-singularity for the phantom case. For 
9 = 2 (a — 0) we have generalized w-singularity. 
The throat radius of the wormhole is 

b(t) 360 



3 - 4ir 2 DgCb {B9- 1 - (B - Ct)^ 1 ) 

es place if and only if ( 
sider the solution ([22 

P (t w ) = -A/3), at 



so the Big Trip takes place if and only if g > 1, i.e., w < — 1 which is the case for the a < 1, as expected. 
4. At last we consider the solution (|22|) with g < 0. Here we have the generalized w-singularity (p(t w ) = 0, 



2V3(\g\ + l)x f 



3A 



1-11-^ 

3/ 



and the Big Trip at 



2^(14+1) _ (l , 



AV3 l\ a:/ 



with 5 = 3\/3/(47rL>6o) > 0. 

Thus, smooth exit from the phantom inflationary phase can still be tentatively recovered by considering a Big 
Trip scenario where the primordial phantom universe would travel in time towards a future universe (filled with, for 
example, usual radiation, see [46|h Such "exit from inflation" is possible in phantom models both with and without 
the future Big Rip singularity. 



Conclusion. In summary, the relatively simple method for constructing phantom energy models without finite-time 
future singularity is developed. The dark energy models without future singularity are attractive from the physical 
viewpoint because the occurrence of finite-time singularity may lead to some inconsistencies. The equivalent descrip- 
tion of the Little Rip cosmology where singularity effectively disappears, via fluid or scalar-tensor theory is presented. 



12 



Phantom models with asymptotically de Sitter evolution are described. It is demonstrated that asymptotically de 
Sitter expansion can be realized in the class of exponential or power-law scalar potentials. Generalization for phantom 
models coupled with dark matter is also discussed. It is interesting to note that disintegration of bound structure 
(on the example of the system of Sun and Earth) in some asymptotically de Sitter phantom universe may occur for 
observationally acceptable choice of parameters. 

We have shown that current data make it essentially impossible to determine whether or not the universe will end 
in a future singularity. The above scalar dark energy models represent natural alternative for ACDM model, which 
also leads to non-singular cosmology. Nevertheless, even for the non-singular asymptotically de Sitter universe, the 
possibility of dramatic rip which may lead to the disappearance of bound structures in the universe remains to be 
possible. 

It is confirmed that phantom models without the Big Rip may lead to wormholes solutions. It is demonstrated 
that the possible Big Trip in the phantom cosmology can happen even if we consider the models with asymptotic de 
Sitter evolution. 

The presence of numbers of free possible parameters (for instance, the choice of equation of state (JSJ)) gives enough 
space for fine-tuning the models which can be useful for fitting with observational data. Hence, the described method 
may be very useful for confronting of theoretical models with coming observational data. 
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